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Water  mass  bio-optical  properties  in  the  Monterey  Bay 
region:  Fluorescence-based  inference  of  shifts 
in  phytoplankton  photophysiology 

J.  K.  Jolliff,1  R.  W.  Gould  Jr.,1  B.  Penta,1  W.  J.  Teague,1  S.  DeRada,1  F.  P.  Chavez,2 
and  R.  A.  Amone1 
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[l]  A  physical  and  bio-optical  field  survey  of  the  Monterey  Bay  area  was  conducted 
during  May-June  2008.  The  combined  bio-optical  and  physical  data  may  be  summarized 
as  a  transition  between  two  end-member  states  during  the  late  spring  to  summer  upwelling 
season:  (1)  the  mesotrophic,  nanoflagellate-dominated,  low-salinity  surface  waters 
(chlorophyll-a  ~  0.5-2  mg  m"3;  S  <  33.4)  of  the  California  Current  and  (2)  the  eutrophic, 
diatomaceous,  higher  salinity  surface  waters  (chlorophyll-a  >  2  mg  m”3;  S  >  33.8) 
of  Monterey  Bay  and  adjacent  continental  shelf  areas.  High-resolution  and  collocated 
spectrophotometric,  fluorometric  and  CTD  data  obtained  from  a  towed  platform  indicated 
low-salinity  subarctic-origin  surface  waters  intruded  into  Monterey  Bay  on  4  June. 

The  dark  in  vivo  fluorometry  (IVF)  phytoplankton  response  normalized  to  particle 
absorption  at  676  nm  (the  apparent  fluorescence  efficiency,  AFE)  was  nearly  fourfold 
larger  in  this  water  mass  type  compared  to  higher  salinity  surface  waters  more  typical  of 
Monterey  Bay.  The  collocated  fluorescence  and  optical  data  were  then  used  to  estimate 
in  situ  irradiance  values  and  determine  apparent  light  saturation  intensities  (rk)  based  on 
the  remarkably  consistent  AFE  water  column  inflection  points.  Ik  values  retrieved  from  the 
low-salinity  surface  waters  were  approximately  half  those  obtained  over  the  continental 
shelf.  An  analysis  of  concomitant  HPLC  data,  in  addition  to  historical  data  for  the  region,  . 
suggest  these  observed  fluorescence  trends  may  be  indicative  of  taxon-specific  variation 
in  photophysiology.  Specifically,  the  subarctic  water  mass-associated  pelagic  nanoflagellate 
group  likely  possesses  a  fundamentally  different  photosynthetic  architecture  than  large 
diatoms  prototypical  of  coastal  upwelling  regimes. 

Citation:  Jolliff,  J.  K„  R.  W.  Gould  Jr.,  B.  Penta,  W.  J.  Teague,  S.  DeRada,  F.  P.  Chavez,  and  R.  A.  Amone  (2012),  Water  mass 
bio-optical  properties  in  the  Monterey  Bay  region:  Fluorescence-based  inference  of  shifts  in  phytoplankton  photophysiology, 

J.  Geophys.  Res .,  //7,  C07019,  doi:  10.1 029/20 1 1JC007568. 


1.  Introduction 

[2]  The  goal  of  the  Naval  Research  Laboratory’s  (NRL) 
Bio-Optical  Studies  of  Predictability  and  Assimilation  for  the 
Coastal  Environment  (BIOSPACE)  program  is  to  forecast 
bio-optical  and  physical  properties  of  the  coastal  zone  on  the 
timescale  of  atmospheric  model  forecasts  (~l-5  days)  using 
coupled  ecosystem-hydrodynamic  numerical  simulations. 
This  goal  requires  an  extensive  interdisciplinary  field  effort 
to  constrain/validate  modeling  activities,  as  well  as  improve 
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our  understanding  of  biological-physical  interactions.  Accord¬ 
ingly,  the  first  BIOSPACE  field  observation  campaign  was 
conducted  during  2008  in  the  Monterey  Bay  region  (Figure  1). 
Following  the  template  of  the  previous  Autonomous  Ocean 
Sampling  Network  (AOSN)  experiments  in  Monterey  Bay 
[Ramp  et  al. ,  2009],  the  BIOSPACE  field  program  sought  a 
comprehensive  collection  of  in  situ  data  from  the  deployment 
of  autonomous  platforms,  moorings,  and  ship  station  sampling 
while  leveraging  available  remote  sensing  information  and 
high-frequency  Doppler  radar  stations. 

[3]  Here  we  present  a  synthesis  of  this  field  effort,  with 
specific  emphasis  upon  relating  biological  and  bio-optical 
properties  with  physical  properties  of  specific  water  mass 
types.  The  concept  of  biological  properties  specific  to  plankton 
being  consistently  delineated  by  independently  described 
physical  water  masses  was  demonstrated  by  Fager  and 
McGowan  [1963]  and  more  recently  applied  to  biological 
studies  of  Monterey  Bay  [Osborn  et  al .,  2007],  Here  we  show 
that  this  is  a  useful  context  for  the  descriptive  analysis  of 
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Figure  1.  (a)  Time  line  of  the  BIOSPACE  ScanFish  (SF)  deployments  (gray)  and  ship  station  sampling 
(red  dots)  with  MBARI  moorings  M2  wind  vectors,  (b)  Map  of  the  SF  and  station  sampling  deployments 
in  the  MB  area.  SF  surveys  on  4  June  arc  shown  in  gray  for  segments  1  and  2.  The  remaining  SF  tracks 
are  in  black.  The  station  sampling  for  the  extraction  of  bottle  chlorophylls,  nutrients,  and  HPLC  data  are 
shown  in  red,  with  a  single  station  where  observed  salinity  was  less  than  33.3  shown  in  blue,  (c)  High- 
frequency  Doppler  radar  estimates  of  surface  currents  in  the  MB  region,  25-h  means  ending  4  June,  0:00 
GMT  are  overlaid  on  the  SF  tracks  (data  available  at  http://www.cocmp.org/). 


our  field  data.  Moreover,  an  understanding  of  the  coupled 
physical  and  bio-optical  features  observed  provides  a  criterion 
for  ecosystcm-hydrodynamic  model  evaluation  that  moves 
beyond  univariate  model  and  data  comparison  statistics  [Jolliff 
et  al. ,  2009;  Stow  et  al. ,  2009]. 

[4]  As  is  typical  of  contemporary  ocean  observational 
campaigns,  the  BIOSPACE  field  study  included  the  deploy¬ 
ment  of  multiple  conventional  chlorophyll-a  fluorometers  on 
a  wide  range  of  platforms  (Slocum  electric  gliders,  towed 
vehicles,  moorings,  autonomous  underwater  vehicles,  CTD 
rosettes,  horizontal  profiling  optics  packages,  underway  ship 
flow-through  systems).  The  2008  field  survey  alone  obtained 
on  the  order  of  ~106  fluorometer  observations.  Ostensibly, 
these  millions  in  vivo  fluorescence  (IVF)  datum  may  serve  as 
a  collective  proxy  for  the  space/time  distribution  of  phyto¬ 
plankton  chlorophyll-tf  concentration  [Lorenzen,  1966]. 
However,  the  IVF  data  collected  appear  to  be  no  better  than 
an  order-of-magnitude  accurate  indicator  of  phytoplankton 
pigment  concentration.  Whereas  these  data  may  be  useful  for 
descriptive  analysis  of  the  presence  or  absence  of  phyto¬ 
plankton  blooms,  more  quantitative  treatment  of  these  data  is 
not  tractable.  This  failing  of  conventional  fluorometry  data  to 
serve  as  a  reliable  proxy  for  phytoplankton  pigment  con¬ 
centration  has  long  been  recognized  [Alpine  and  Cloem , 
1985;  Cullen ,  1982;  Falkowski  and  Kolber ,  1995;  Loft  us 
and  Seliger,  1975;  Strickland,  1968]. 

[5]  Nevertheless,  fluorescence  is  a  specific  physiological 
response  of  the  photosynthetic  apparatus;  effectively  a 


pathway  of  activated  chlorophyll  deexcitation  in  competition 
with  the  photochemical  pathway  and  energy  dissipation  as 
heat  [Krause  and  Weis ,  1991].  An  alternative  usage  of  IVF 
data  is  to  place  the  observed  fluorescence  intensity  response  in 
the  context  of  the  overall  physiological  energy  transformation 
within  photosystem  II  (PSI1)  of  the  phytoplankton  photosyn¬ 
thetic  apparatus  [see  Huot  and  Babin ,  2010].  Fluorescence 
properties  observed  in  situ  are  thus  indicative  of  photo- 
physiological  properties  of  natural  phytoplankton  populations 
[Falkowski  and  Kolber ,  1995;  Kiefer  and  Reynolds ,  1992]. 
The  challenge  presented  by  these  IVF  data  is  to  carefully  make 
inferences  about  such  photophysiological  properties  via  anal¬ 
ysis  of  the  additional  bio-optical  and  environmental  variables 
observed  [< Cullen ,  1982;  Kiefer ,  1973;  Sackmann  etal .,  2008]. 
Herein  IVF  data  are  examined  where  concomitant  spectral 
absorption,  scattering,  and  backscattcring  coefficients  were 
observed  as  well  as  the  physical  properties  of  seawater.  Inde¬ 
pendent  estimates  of  phytoplankton  pigment  concentration, 
particle  absorption,  and  light  attenuation  provided  an  appro¬ 
priate  context  to  make  inferences  from  these  IVF  data  that  may 
improve  understanding  of  integrated  biological-physical  pro¬ 
cesses  in  this  region. 

2.  Methods 

2.1.  Study  Site 

[6]  The  Monterey  Bay  (hereinafter  MB)  region  is  part  of 
the  California  Current  System  (CCS),  and  it  is  located  within 
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an  extended  regional  transition  zone  between  the  subarctic 
Pacific  and  the  North  Pacific  subtropical  gyre  [Roden, 
1991].  The  subarctic  front,  approximately  the  33.8  isoha- 
linc  [Tully  and  Barber ,  I960],  lies  along  ~40-43°N  latitude 
and  curves  southward  west  of  150°W  to  form  the  western 
boundary  of  the  CCS  [Lynn,  1986].  Thus,  water  mass  char¬ 
acteristics  of  the  transition  zone  are  influenced  by  the  dif¬ 
ferent  temperature  and  salinity  (T/S)  properties  of  the 
subarctic  and  the  subtropics.  In  the  MB  region,  apparent  T/S 
signatures  of  both  water  mass  types  may  be  found  [  Warn - 
Vamas,  2007].  The  relative  predominance  of  these  con¬ 
trasting  water  mass  characteristics  at  any  given  time  and 
location  is  brought  about  by  a  complex  interplay  between 
the  meandering  and  shoreward  translation  of  the  equator- 
ward  flowing  California  Current  [Checkley  and  Barth, 
2009],  the  prevailing  mode  of  the  wind  stress  (upwclling/ 
downwelling-favorable)  [Storlazzi  et  al,  2003],  mesoscale 
eddies  and  their  interaction  with  both  the  California  Current 
and  upwelling  jets  [Paduan  and  Rosenfeld,  1 996],  and  the 
seasonally  recurring  surface  expression  of  poleward  flows 
along  the  California  coast  [Chelton  et  al.,  1988]. 

[7]  Herein  an  attempt  is  made  to  relate  these  physical 
water  mass  characteristics  to  bio-optical  signatures.  The 
pelagic/low-salinity  subarctic  Pacific  is  regarded  as  an  iron- 
limited  system  [Boyd  et  al.,  2004],  and  in  contrast  to  similar 
latitudes  in  the  North  Atlantic,  there  is  no  distinct  spring 
bloom.  Subarctic  Pacific  surface  chlorophyll-0  concentration 
(hereinafter  chi)  values  there  arc  consistently  on  the  order  of 
~0.5-2  mg  m  3  [Miller  et  al,  1991].  Upwelling  along  the 
central  California  coast,  however,  results  in  surface  chi 
values  an  order  of  magnitude  higher  [Collins  et  al,  2003]. 
Within  MB  (an  open  embayment  spanning  approximately 
40  km  north-to-south),  a  highly  diverse  and  variable  phyto¬ 
plankton  assemblage  is  observed  [Rines  et  al,  2010]  and  chi 
values  range  to  >500  mg  m  3  [Lee  et  al,  2007;  Ryan  et  al, 
2009]. 

[8]  With  respect  to  temperature  and  salinity,  California 
Current  water  tends  to  be  cold  and  fresh  due  to  a  subarctic 
origin,  whereas  MB  surface  waters  and  surface  waters  of 
more  southerly  origin  arc  warmer  and  more  saline,  reflecting 
near-equatorial  sources  [Breaker  and  Broenkow,  1994], 
Accordingly,  Spiciness  (n)  is  a  useful  state  variable  for  this 
analysis  [Flament,  2002].  n  isoplcths  arc  approximately 
orthogonal  to  those  of  density  in  T/S  space  [Fofonoff,  1985] 
such  that  on  an  isopycnal  surface  some  water  masses  may  be 
“spicier”  (warmer/saltier)  than  others.  This  analysis  is 
largely  restricted  to  the  upper  ~80  m  during  a  seasonal 
period  of  upwelling  favorable  wind  stress  along  the  coast 
(Figure  la);  diapyncal  mixing  of  subarctic-origin  California 
Current  surface  waters  with  recently  upwelled,  higher 
salinity  waters  near  the  coast  appears  to  generally  follow 
constant  n  isopleths.  Large  n  deviations  from  the  prevailing 
n  isoplcths  were  interpreted  as  indicative  of  the  addition 
(positive-n)  or  removal  (negative-FI)  of  thermal  energy  due 
to  air/sea  exchange  processes. 

2.2.  Bio-Optical  Survey  and  Analysis  Methods 

[9]  As  part  of  NRL’s  BIOSPACE  project,  an  extensive 
field  campaign  was  conducted  in  the  MB  area  from  30  May 
to  19  June  2008  based  primarily  on  board  the  R/V  Point  Sur 
(Moss  Landing  Marine  Laboratories).  Observations  included 
ship  station  CTD  casts  (Figure  1)  with  multiple  analyses 


performed  on  water  samples  from  the  Niskin  bottle  rosette. 
Macronutrient  and  chlorophyll  assays  were  performed  using 
methods  detailed  in  Pennington  and  Chavez  [2000].  High 
Performance  Liquid-Chromatography  (HPLC)  analyses  were 
performed  on  near-surface  water  samples  filtered  through 
GF/F  filters,  stored  in  liquid  nitrogen  and  shipped  to  the  Horn 
Point  Laboratory  (University  of  Maryland)  for  analysis  fol¬ 
lowing  standard  NASA  HPLC  protocols  [Hooker  et  al, 
2005].  Aliquots  of  GF/F-filtered  and  liquid  nitrogen-stored 
sample  water  were  also  obtained  for  spectrophotometric 
determination  of  absorption  coefficients  (Analytical  Spectral 
Devices  FicldSpcc  spcctroradiomcter  with  a  spectral  range 
of  ~338  nm-~1069  nm  at  ~1.4  nm  resolution).  Total 
absorption  coefficients  were  partitioned  into  phytoplankton, 
particulate  detritus,  and  dissolved  contributions  following 
methods  described  in  Mitchell  et  al  [2000]. 

[10]  This  analysis  focuses  primarily  on  the  observations 
from  the  ScanFish  MK  II  (SF)  platform  (Figure  1).  The  SF  is 
a  remotely  operated  and  towed  vehicle  that  was  piloted  in  the 
water  column  from  just  below  the  surface  to  about  80  m 
depth.  A  multiple  sensor  suite  was  mounted  on  the  SF  plat¬ 
form:  AC-9  spectrophotometer  for  the  measurement  of  mul- 
tispectral  absorption  and  scattering  coefficients;  WET  Labs 
Environmental  Characterization  Optics  (ECO)-BB3  instru¬ 
ment  used  here  for  the  retrieval  of  particle  backscattering 
coefficients,  Sea-Bird  Electronics  CTD,  and  a  WET  Labs 
ECO  chlorophyll-0  fluorometcr.  Here  we  isolate  2170  quasi- 
vertical  profiles  from  both  the  descent  and  ascent  of  the 
towed  platform,  with  a  mean  spatial  resolution  of  ~2  obser¬ 
vations  per  vertical  meter.  The  AC-9  raw  data  were  corrected 
for  temperature,  salinity,  and  scattering  effects  [Gould  et  al, 
1999;  Pegau  et  al,  1997;  Zaneveld  et  al,  1994]. 

[11]  The  chlorophyll  fluorescence  sensor  was  factory  cal¬ 
ibrated  to  yield  apparent  “chlorophyll”  values  using  a  linear 
operator  and  these  data  were  not  corrected  for  additional 
quenching  effects  [Cullen,  2008;  Sackmann  et  al,  2008]. 
Instead,  the  fluorometer  chlorophyll  concentration  values 
were  reassigned  as  arbitrary  fluorescence  units  (AFU),  i.e., 
these  data  were  interpreted  as  an  indicator  of  the  fluorescence 
intensity  response  without  additional  presumptions  regarding 
pigment  concentration  or  phytoplankton  biomass.  Similar 
treatment  of  conventional  “chlorophyll-0”  fluorometer  data 
as  “relative  fluorescence  units”  may  be  found  in  Ryan  et  al 
[2005]  and  Yu  et  al  [2002]. 

[12]  Data  from  the  SF-mounted  AC-9  spectrophotometer 
were  used  as  an  alternative  proxy  for  phytoplankton  pigment 
concentration.  An  AC-9-based  particle  absorption  at  676  nm 
(0p-676)  was  calculated  by  subtracting  the  literature  pure  water 
value  [Pope  and  Fry,  1997]  from  the  total  absorption  coeffi¬ 
cient  observations  and  presuming  the  absorption  contribution 
from  non-living  organic  matter  was  negligible  at  longer 
wavelengths.  The  observed  relationship  between  our  spectro¬ 
photometric  filter  pad  measurements  of  particle  absorption 
against  the  bottle  chi  measurements  provided  an  empirical  fit 
that  was  applied  to  the  AC-9  data  set  (Figure  2).  The  robust 
correspondence  between  log-transformed  0p-676  and  bottle 
chi  determinations  (r2  =  0.88,  n  =  160)  was  consistent  with 
previous  studies  [Babin  et  al,  2003;  Bricaudet  al,  1995]  and 
provided  a  means  to  transform  the  SF-mounted  AC-9  data  into 
a  new  proxy  for  phytoplankton  pigment  concentration,  which 
is  referred  to  hereinafter  as  chl-0p  (mg  m“3). 
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Figure  2.  Bottle  chlorophyll  measurements  from  ship  sta¬ 
tion  sampling  versus  speetrophotometrie  determination  of 
particle  absorption  coefficients  (676  nm). 


eaeh  individual  sample  [Suggett  et  al .,  2009].  Here,  the 
analysis  compensates  for  an  incomplete  fluorescence  diag¬ 
nostic  on  each  datum  by  examining  trends  for  a  very  large 
aggregate  “weak  probe”  1VF  data  set  (~2  x  105)  in  the 
context  of  concomitantly  collected  optical  and  physical  data. 
Particularly  useful  was  the  estimate  of  the  in  situ  irradiance 
field  integrated  over  the  Photosynthetieally  Active  Radiation 
spectral  range  (400-700  nm).  These  irradiance  estimates 
were  calculated  by  combining  the  ship  meteorological  data 
with  the  SF  profile  Inherent  Optical  Property  (IOP)  data. 
Details  for  this  procedure  are  provided  in  section  Al. 

[15]  The  ancillary  “historical  data”  used  in  this  analysis 
refers  to  the  bimonthly  California  Cooperative  Fisheries 
Investigations  (CalCOFI)  Line  67  station  surveys  conducted 
between  1989  and  2007  and  biweekly  hydrographic  and 
biological  surveys  conducted  near  MBARI  permanent 
mooring  Ml  (36°44.8,N,  122°1.3W)  [Chavez  et  al ,  2002; 
Pennington  and  Chavez ,  2000].  These  data  were  collected 
and  archived  by  the  Monterey  Bay  Aquarium  Research 
Institute  (MBARI)  Ocean  Biology  Group,  and  post-1996 
biweekly  cruises  were  conducted  as  part  of  the  Studies  of 
Ecological  and  Chemical  Responses  to  Environmental 
Trends  (SECRET)  program  [Chavez  et  al ,  2002].  Further 
information  on  data  retrieved  from  MBARI  permanent 
moorings  and  previous  mieronutrient  surveys  of  the  MB  area 
are  provided  in  section  A2. 


[13]  When  the  AC-9-bascd  proxies  for  phytoplankton 
pigment  concentration  were  compared  to  the  collocated  SF 
fluorescence  data,  patterns  emerged  that  were  broadly  con¬ 
sistent  with  what  is  known  about  phytoplankton  fluores¬ 
cence  physiology.  Since  it  is  not  customary  to  analyze 
conventional  chlorophyll  fluorometer  data  in  this  manner, 
some  additional  variables  were  defined  to  aid  in  the  analysis. 
The  Apparent  Fluorescence  Efficiency  (AFE)  was  defined  as 
the  in  vivo  fluorescence  response  (F,  Ex/Em  470/695  nm, 
quantified  as  AFU)  normalized  by  the  AC-9  chlorophyll 
estimate  (ehl-ap) 


AFEC  = 


F(AFU) 

chl-ap 


(1) 


The  subscript  “e”  indicates  that  the  AFE  value  is  determined 
from  the  chlorophyll  estimate  based  on  the  AC-9  absorption 
value.  Similarly,  the  AFE  values  were  also  examined  with 
respect  to  the  original  op-676  nm  coefficient 


AFEp  = 


F(AFU) 
ap-676  nm’ 


(2) 


[14]  In  both  calculations  the  intent  is  to  estimate  the  fluo¬ 
rescence  efficiency  of  the  phytoplankton  population.  The  word 
“apparent”  is  used  throughout  to  emphasize  that  conven¬ 
tional  ehl  fluoromteiy  (i.e.,  excitation  with  a  weak  flash  of 
light  that  does  not  significantly  perturb  the  electron  transport 
system  of  the  photosynthetic  units)  is  quite  distinct  from 
“active”  or  “variable”  fluorometry  techniques  such  as  Fast 
Repetition  Rate  Fluorometry  (FRRF)  wherein  repetitive  or 
saturating  flashes  of  light  arc  used  to  fully  reduce  photo¬ 
system  II  (PSII)  to  obtain  the  maximum  fluorescence 
response  and  other  derivative  fluorescence  parameters  on 


3.  Results 

3.1.  Bio-Optical  Properties  and  Water  Mass  Types 

[16]  The  SF  platform  tows  were  first  conducted  during 
4  June  along  a  transect  beginning  within  MB  (Segment  1; 
Figure  lb)  and  resuming  outside  of  MB  (Segment  2)  during  a 
brief  pause  in  the  prevailing  northwesterly  winds  (Figure  1  a). 
The  potential  shoreward  intrusion  of  subarctic-origin  water 
was  indicated  by  SF  surface  salinities  below  33.4  (Segment  2; 
Figure  3);  whereas  within  MB  the  SF-detected  salinities 
remained  above  33.8  (Segment  1;  Figure  3).  Low  temperature 
(~9.0°C)  and  comparatively  high  salinity  (~34.0)  observa¬ 
tions  made  within  MB  (Segment  1)  were  likely  indicative  of 
recently  upwelled  North  Pacific  Deep  Water  (NPDW).  Higher 
temperatures  for  salinities  >33.8  (a  positivc-Fl  deviation) 
suggest  subsequent  warming  of  NPDW  within  MB.  Warmer 
waters  were  coincident  with  higher  total  440  nm  absorption 
coefficients  in  both  transects  (Figure  3a),  although  MB  surface 
water  absorption  coefficients  were  an  order-of-magnitude 
higher. 

[17]  The  4  June  T/S  observations  on  Segment  2  appear  to 
have  captured  atypical  conditions  with  respect  to  historical 
observations  in  this  area.  Only  13.4%  of  the  historical  CTD 
values  (n  =  55,015)  recorded  in  the  upper  200  m  meet  a 
negative-n  criterion  over  the  25  <cr0<  26  density  interval,  in 
contrast  to  60%  of  the  Segment  2  SF  observations.  Coinci¬ 
dent  Mooring  M2  observations  confirm  the  presence  of  this 
negative-n  water  mass  in  the  area  (Figure  4).  A  review  of 
these  M2  T/S  data  (1999-2008)  revealed  that  the  frequency 
of  appearance  for  this  water  mass  type  ranged  from  ~1— 30% 
of  M2  observations  during  the  May  through  June  period  of 
each  year. 

[is]  This  low-salinity,  negative-n  water  mass  type 
appeared  to  be  a  cold  haloeline  form  of  subarctic  upper 
water  (SUW)  and  was  likely  adveeted  into  the  MB  area  via 
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Figure  3.  T/S  plot  of  4  June  SF  surveys,  Segments  1  (S  > 
33.8)  and  2  (S  <  33.8).  Color  contours  are  the  AC-9  total 
absorption  coefficients  (440  nm)  at  each  point.  Density  and 
n  (dashed  line)  contours  are  overlaid  in  T/S  space. 

the  meandering  of  the  offshore  California  Current  [Wam- 
Varnas ,  2007].  Indeed,  CODAR  surface  current  estimates 
indicated  the  presence  of  an  anticyclonic  eddy  centered  just 
west  of  Mooring  M2  at  the  time  of  the  SF  Segment  2 
observations  (Figure  lc).  Negative  Spiciness  is  interpreted 
here  as  an  indication  of  a  comparative  thermal  energy  deficit 
over  the  observed  salinity  range.  The  sign  of  the  Spiciness 
state  variable  is  an  arbitrary  convention;  however,  it  conve¬ 
niently  contrasts  with  more  prototypical  positive  halocline  n. 
Previous  observations  of  atypically  negative-n  upper  halo¬ 
cline  water  in  the  CCS  have  been  attributed  to  anomalous 
surface  cooling  in  the  subarctic  source  water  region  followed 
by  advection  farther  south  [Freeland  et  al .,  2003]. 

[19]  The  remainder  of  the  SF  survey  T/S  properties 
generally  fell  into  one  of  three  T/S  domains:  (1)  higher 
salinity  MB  surface  water  that  appears  to  result  from  the 
mixing  and  warming  of  recently  upwellcd  NPDW;  (2)  low- 
salinity,  positive-n  surface  water  from  the  offshore  transect 
(OT,  Figure  lc);  and  (3)  transitional  T/S  surface  water 
properties  between  these  end  points.  The  33.2  to  33.8  halo¬ 
cline  waters  for  the  remainder  of  the  data  were  confined  to  the 
0.0  to  0.5  n  isopleths  (Figure  4).  These  data  suggest  that 
negative-n  SUW  is  distinct  from  typical  CCS  surface  waters 
in  this  region.  One  may  interpret  the  more  typical  warmer 
halocline  waters  as  indicative  of  diapycnal  mixing  between 
upwelling  of  colder,  higher  salinity  waters  along  the  shelf  and 
higher  temperature  waters  within  the  coastal  to  California 
Current  “transition  zone”  [Collins  et  al .,  2003]  where  iso- 
pycnal  surfaces  slope  upwards  toward  the  margin  and  the 


33.4  surface  salinity  isopleth  delineates  coastal  waters  from 
the  offshore  California  Current. 

[20]  This  analysis  of  water  mass  physical  properties  sug¬ 
gests  that  the  phytoplankton  population  of  Segment  2  (nega¬ 
tive-n;  S  <  33.4)  is  likely  indicative  of  a  truly  pelagic  and 
subarctic  character  without  diapycnal  mixing  with  shoreward 
water  masses,  and  by  inference,  coastal  phytoplankton  popu¬ 
lations  found  there.  Moreover,  Segment  1  and  2  fluorescence 
data  were  recorded  at  night,  thus  sunlight-stimulated  fluores¬ 
cence  quenching  was  unlikely  to  be  a  contaminant  of  these 
data.  Indeed,  the  fluorescence  signal  per  unit  absorption  from 
Segment  2  was  over  threefold  higher  than  Segment  1 
(Figure  5).  Although  the  total  absorption  coefficient  ranges  for 
Segment  1  (representative  of  MB  surface  waters)  ranged  to 
much  larger  values,  the  apparent  increase  in  AFEp  for  Seg¬ 
ment  2  was  evident  even  where  the  particle  absorption  values 
were  comparable  (e.g.,  0.02  m“!  ap  676;  Figure  5a).  Conver¬ 
sion  of  these  ap-676  values  to  chl-crp  revealed  the  same  trend 
with  respect  to  the  fluorescence  response;  however,  maximum 
Segment  2  surface  chl-ap  values  were  in  the  ~2-3  mg  m“3 
range,  whereas  Segment  1  surface  chl-ap  values  ranged 
upwards  to  ~15  mg  m“3.  This  suggests  that  the  AFE  differ¬ 
ences  may  be  concomitant  to  a  transition  in  trophic  state 
between  the  respective  water  masses,  i.e.,  mesotrophy  (~2  mg 
m“3  surface  chi)  to  eutrophy  (>2  mg  m”3  surface  chi 
[cf.  Bricaud  et  al.,  2004]). 

[21]  If  indeed  dark  AFE  is  a  distinguishable  bio-optical 
signature  of  phytoplankton  populations  that  are  endemic  to 
specific  water  mass  types,  then  it  would  be  expected  that 
Segments  1  and  2  may  serve  as  end  points  with  dark  AFE’s 
at  intermediate  salinities  ranging  quasi-conservatively 
between  them.  This  expected  trend  was  observed  (Figure  6). 
These  AFEC  observations  for  dark  (<10  /.iE  m-2  s-1  incident 


33.0  33.2  33.4  33.6  33.8  34.0  34.2 

Salinity 


Figure  4.  All  SF  data  are  shown  in  black  with  the  4  June 
SF  Segment  2  shown  in  gray.  The  mooring  M2  observations 
(red)  confirm  the  negative-n  water  mass  presence  also 
detected  by  the  SF-mounted  CTD. 
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Figure  5.  (a)  Fluorescence  (F;  AFU)  versus  AC-9  esti¬ 
mated  particle  absorption  at  676  nm  for  Segment  1  (black) 
and  Segment  2  (gray):  the  linear  regression  line  is  shown 
in  red  and  the  slope  is  indicated  as  AFEp.  (b)  As  in  Figure  5a 
but  the  chl-arp  values  are  used  instead  and  the  AFEC  slope 
value  estimates  are  given  for  the  respective  transects. 


scalar  PAR)  and  near-surfacc  waters  (<20  m  depth)  appeared 
to  corroborate  the  notion  that  certain  biological  properties 
of  plankton  may  be  delineated  by  independently  described 
physical  properties  of  seawater  [Fager  and  McGowan,  1 963]. 

3.2.  Apparent  Fluorescence  Quenching 

[22]  Exceptions  to  these  AFE  patterns  occurred  where  near¬ 
surface  fluorescence  observations  were  made  during  daylight 
hours.  For  example,  the  fluorescence  relationships  to  arp-676 
for  the  night  phase  of  the  outer  transect  (OT,  Figure  lc)  were 
approximately  log  linear  (Figure  7a).  This  relationship  was 
used  as  an  initial  proxy  for  the  expected  dark  fluorescence 
response,  i.e.,  it  was  assumed  that  the  PSII  reaction  centers 
were  open  (represented  by  the  term:  F*o)  during  this  time. 
Quenching  of  fluorescence  was  inferred  by  deviation  from  this 
tfp-676  -to-F*o  log  linear  relationship  during  daylight  hours. 
The  deviation  from  dark  values  was  positive  toward  an 
apparent  peak  in  fluorescence  and  then  a  significant  negative 
deviation  occurred  (Figure  7b). 

[23]  Since  fluorescence  quenching  was  suggested  by 
deviations  from  F*o,  the  estimate  of  F*o  was  subtracted 
from  the  observed  fluorescence  intensity  (F;  AFU)  to  cal¬ 
culate  an  apparent  variable  fluorescence  (F*v) 

F*v  =  F  —  F*o,  (3) 

where 

F*o  =  io(,oglo(<,p_676)X+Y)  (4) 

and  the  log  linear  regression  coefficients  (X,  Y)  were 
determined  from  the  nearest  set  of  available  dark  fluores¬ 
cence  and  tfp-676  observations. 

[24]  The  F*v  profiles  resulting  from  these  calculations 
mimicked  fluorescence  induction  curves:  F*v  increased  from 
near  zero  at  depth  to  a  local  water  column  peak  as  PS11 
reaction  centers  were  presumably  closed  and  fluorescence 
efficiency  increased  (photochemical  quenching  variability; 
PQ,  Figure  7c).  Ascending  higher  into  the  water  column,  the 
inflection  back  toward  lower  F*v  values  was  then  likely  due 
to  the  dominance  of  nonphotochemical  quenching  (NPQ), 
i.e.,  high-light  intensities  resulted  in  potential  oxidative 
stress  or  damage. 

[25]  It  is  important  to  distinguish  the  apparent  variable 
fluorescence  (F*v)  peak  observed  here  from  the  fluorescence 
maximum  as  observed  with  variable  fluorescence  techniques 
[Falkowski  and  Kolber ,  1995;  Kolber  and  Falkowski ,  1993]. 
The  true  florescence  maximum  occurs  when  all  the  primary 
quinone-type  electron  acceptors  (QA)  of  PSII  are  fully 
reduced  [Krause  and  Weis ,  1991],  and  this  may  not  have 
occurred  in  the  apparent  maximum  of  these  F*v  profiles. 
Nevertheless,  it  is  presumed  that  the  F*v  inflection  point 
under  ambient  irradiance  indicates  the  effective  saturation  of 
PSII  reaction  centers;  thus  the  F*v  inflection  points  were 
interpreted  as  apparent  irradiance  saturation  values  (1'k).  These 
values  for  the  outer  transect  profiles  were  remarkably  consis¬ 
tent  at  ~67  /iE  m  2  s~ 1  (Figure  8a).  The  zero  crossing,  where 
the  observed  fluorescence  was  below  the  expected  dark  fluo¬ 
rescence,  consistently  occurred  at  ~405  fiE  m'2  s_1.  Specu¬ 
latively,  negative  F*v  values  may  have  been  indicative  of 
true  photoinhibitory  damage  to  PSII  reaction  centers,  i.e., 
reaction  centers  damaged  from  oxidative  stress  were  no  longer 
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Figure  6.  Estimates  of  Apparent  Fluorescence  Efficiency  (AFEC)  versus  salinity  for  all  dark  and  near  sur¬ 
face  (<20  m)  SF  observations:  Gray  asterisks  correspond  to  the  4  June  Segment  2  observations  (top  left) 
and  the  Segment  1  observations  (bottom  right),  all  other  observations  are  black.  The  linear  regression 
(AFEC  versus  salinity)  trend  line  is  shown  (dashed  line). 


functional  [Maxwell  and  Johnson,  2000].  This  is  in  contrast  to 
eneigy-dependent  nonphotochemical  quenching  that  involves 
changes  in  the  PS II  absorption  cross-section  and  thermal  dis¬ 
sipation  of  energy  but  is  not  associated  with  structural  damage 
to  the  photosynthctic  apparatus  [Kiefer  and  Reynolds ,  1992; 
Maxwell  and  Johnson ,  2000].  AFEp  versus  in  situ  irradiance 
patterns  confirm  the  identical  apparent  irradiance  saturation 
peak  at  67  pE  m“2  s_l  (Figure  8b). 

[26]  Whereas  the  fluorescence  versus  irradiance  (F  versus  7) 
curves  are  not  functionally  identical  to  photosynthesis  versus 
irradiance  curves  (P  versus  I)  [Kolber  and  Falkowski ,  1993; 
Stegmann  et  al.,  1 992],  the  apparent  F*v  inflection  points  here 
are  likely  to  occur  at  an  in  situ  irradiance  close  to  the  saturation 
of  photosynthctic  activity  [Kiefer  and  Reynolds ,  1992].  These 
inferred  Ik  values  were  low  compared  to  literature  (P  versus  I 
curve)  Ik  values  [Jassby  and  Platt ,  1976]  and  would  be 
indicative  of  dark-acclimated  microflora  [Chan,  1978],  or  Ik 
values  more  typical  of  high-latitude  phytoplankton  [Gilstad 
and  Sakshaug,  1990].  The  genuine  Ik  values  may  be  differ¬ 
ent  than  these  estimated  quantities,  but  it  was  apparent  that 
phytoplankton  endemic  to  low-salinity/subarctic-origin  sur¬ 
face  water  mass  types  were  consistently  responding  to  oxida¬ 
tive  stress  at  comparatively  low  irradiances  throughout  the 
solar  period  of  the  offshore  SF  transect  (Figures  9a  and  9b). 
The  surface  salinity  was  below  33.4  and  surface  chi  was  on  the 
order  of  ~2  mg  m  3  (Figure  9c),  suggesting  that  this  “off¬ 
shore”  transect  was  representative  of  California  Current 
mcsotrophic  conditions. 

[27]  These  F*v  versus  irradiance  calculations  were 
repeated  for  a  case  of  apparent  nonphotochemical  quenching 
in  surface  waters  much  closer  to  the  coast  and  over  the  con¬ 
tinental  shelf  northwest  of  MB  (ST  in  Figure  1  c).  The  mean 
surface  salinity  for  this  transect  was  34.0;  mean  surface  chi 
(chl-ap)  was  3.3  mg  m  ~3  and  ranged  as  high  as  6.9  mg  m~3. 


Here  again  the  apparent  quenching  patterns  (as  a  function 
of  the  estimated  in  situ  irradiance)  were  consistent.  How¬ 
ever,  the  collective  “shelf  water”  patterns  found  here  were 
distinct  from  those  observed  offshore.  The  inflection  point 
(Ik)  was  observed  at  irradiant  intensities  more  than  double 
the  offshore  values:  ~  147  pE  m  2  s-1.  Further,  the  F*v 
values  did  not  fall  below  zero  over  the  same  high  irradiance 
levels  (>400  pE  m“2  s-1;  Figure  10a).  Peak  AFEp  values 
occurred  at  twice  the  in  situ  irradiance  as  the  offshore  tran¬ 
sect,  but  the  peak  AFEp  values  were  half  as  high.  AFEp 
subsurface  peaks  were  also  detected  on  some  SF  profiles 
obtained  within  MB  (data  not  shown)  at  a  similar  in  situ 
irradiance  value  (~1 52  pE  m“2  s"*1).  This  suggests  there  is 
some  inherent  distinction  between  neritic  phytoplankton 
populations  and  those  farther  offshore  that  is  manifest  in 
these  apparent  irradiance  saturation  values. 

3.3.  Phytoplankton  Species  Composition  and  Dissolved 
Nutrients 

[28]  Bulk  contrasts  in  the  phytoplankton  community  fluo¬ 
rescence  responses  may  be  due  to  differences  in  phytoplank¬ 
ton  community  composition  and  associated  variations  in 
pigment  content  [Alpine  and  Cloem,  1985;  Beutler  et  al., 
2003;  Falkowski  and  Kolber ,  1995]  as  well  as  nutritional  sta¬ 
tus  [Behrenfeld  et  al,  2006;  Hou  et  al,  2007].  Both  of  these 
conditions,  in  tum,  may  impact  light  acclimation  responses 
and  associated  changes  in  apparent  fluorescence  properties. 
Strictly  taxon-specific  apparent  fluorescence  property  differ¬ 
ences  would  be  consistent  with  the  notion  that  subarctic-origin 
surface  water  masses  contain  pelagic  pico-  and  nanophyto- 
plankon  populations  (<20  pm  cell  diameter)  whereas  MB 
waters  contain  neritic  populations  composed  predominantly  of 
diatomaceous  microphytoplankton  (>20  cell  diameter) 
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[Chavez  et  al. ,  1991;  Collins  et  al. ,  2003;  Wilkerson  et  al , 
2000]. 

[29]  HPLC  data  collected  during  station  sampling  allow 
for  an  evaluation  of  potential  fluorcsccnce-based  inferences 
of  phytoplankton  community  shifts  and  associated  differ¬ 
ences  in  pigment  content.  The  mean  pigment-based  micro¬ 
phytoplankton  proportion  factor  [Hooker  et  al ,  2005]  was 
90%  (n  =  58)  during  the  station  survey,  with  one  notable 
outlier  at  30%  (Table  1).  This  reflects  a  station  sampling 
bias  confined  closer  to  the  coast  and  the  interior  of  MB 
(Figure  lb).  Nonetheless,  the  respective  outlier  bottle  salinity 
was  33.28,  compared  to  S  >  33.70  for  the  remaining  HPLC 
samples.  Significant  departures  of  this  outlier  sample  from 
the  expected  HPLC  variable  distributions  (p  <  0.01)  occurred 
for  a  greater  proportion  of  chlorophyll-6  and  photoprotcctive 
carotenoids  (Table  1).  The  sample  was  also  significantly  ele¬ 
vated  in  19'-butaloyloxyfucoxanthin,  19-hexanoyloxyfucox- 
anthin,  and  lutein;  generally  indicative  of  greater  pclagophyte, 
prymnesiophyte,  and  chlorophyte  abundance,  respectively 
[Riegman  and  Kraay ,  2001].  This  group  is  generally  <20  /*m 
cell  diameter  and  may  be  referenced  as  the  nanoflagellate 
group  [Jeffrey  and  Vesk ,  1997].  Pigment  markers  for  this 
nanoflagellate  group  are  a  distinctive  feature  of  HPLC  data 
obtained  in  the  subarctic  North  Pacific  [Fujiki  et  al ,  2009]. 

[30]  Similar  trends  were  evidenced  by  alternative  taxo¬ 
nomic  identification  methods  in  the  historical  data.  Cell 
counts  collected  over  1989-2006  were  sorted  by  total  auto¬ 
trophic  biomass  and  the  mean  percentage  contributions  from 
picophytoplankton  and  diatoms  were  calculated  (Figure  1 1  a). 
Diatoms  crossed  50%  of  the  estimated  total  chi  threshold  at 
approximately  ~2  mg  m-3  chi.  Thus  the  threshold  between 
mcsotrophic  and  cutrophic  surface  waters,  established  herein 
at  ~2  mg  m~3  total  chi,  was  also  an  approximate  transition 
point  into  dominance  of  diatomaccous  flora.  These  percent¬ 
age  estimates  were  calculated  by  converting  cell  volume 
carbon  estimates  to  chlorophyll-^  values  via  a  constant 
conversion  factor  of  47.2,  which  is  the  mean  carbon-to- 
chlorophyll-tf  ratio  when  cell  carbon  estimates  were  com¬ 
pared  to  bottle  chlorophyll-^  assays  in  the  historical  data 
(Figure  lib).  Collectively,  these  data  are  consistent  with 
previous  work  demonstrating  a  general  trend  toward  micro- 
phytoplankton  dominance  at  cutrophic  (>2  mg  m  3)  chi 
values  [Bricaud  et  al ,  2004;  Li,  2002],  which  tend  to  consist 
of  large  diatoms  in  coastal  upwclling  systems  [Walsh,  1988]. 

[31]  In  the  SF  data  set,  the  low-salinity  phytoplankton 
communities  (S  <  33.4)  reached  maximum  chl-ap  values  of 
~2-3  mg  m  3;  higher  chl-ap  values  corresponded  to  higher 
salinities  and  a  shift  in  the  IVF  responses  toward  lower 
AFE’s.  Given  the  hypothesis  that  phytoplankton  community 
species  composition,  cell  size,  and  pigment  composition 


Figure  7.  Log  plot  of  fluorescence  (F,  AFU)  versus  particle 
absorption  at  676  nm  (a)  during  the  night  hours  (01:52  -04:51 
local  time)  of  the  outer  transect  and  (b)  during  daylight  hours 
(13:21-14:43  local  time).  The  expected  “dark”  fluorescence 
response  regression  line  as  a  function  of  ap-676  nm  is  shown 
in  Figures  7a  and  7b.  (c)  Daylight  hour  profiles  of  fluorescence 
from  the  outer  transect  (13:21—14:43  local  time).  The  apparent 
variable  fluorescence  (F*v)  calculation  is  described  in  the  text 
(PQ)  indicates  presumed  photochemical  quenching  and  (NPQ) 
indicates  presumed  nonphotochemical  quenching. 
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Figure  8.  (a)  Apparent  variable  fluorescence  (F*v)  versus  irradiance  for  the  outer  SF  transect  survey. 
The  trend  line  (red)  is  the  smoothed  data,  dashed  lines  indicate  the  inflection  point  and  zero  crossing 
of  the  trend  line,  (b)  The  AFEp  for  the  same  set  of  points;  the  dashed  line  indicates  the  inflection  point 
for  the  smoothed  data. 


generally  trend  in  broadly  consistent  ways  as  a  function 
of  the  total  autotrophic  biomass  in  this  system,  nutrient 
concentrations  were  also  examined  as  an  additional  con¬ 
straint  on  the  transition  from  pelagic  mesotrophy  to  coastal 
eutrophy.  As  observed  elsewhere  [Kamykowski  etal,  2002], 
nitrate  versus  temperature  plots  for  the  upper  ocean  (<200  m) 
conformed  to  a  linear  trend  in  our  data  as  well  as  histor¬ 
ical  data  for  the  MB  region  (Figure  12).  The  nitrate  versus 
temperature  zero  ordinate  intercept,  the  nitrate  depletion 
temperature  (NDT),  was  14.2°C  during  the  2008  BIO¬ 
SPACE  survey  (n  =  638;  r2  =  0.89)  and  14.4°C  for  the  his¬ 
torical  data  (n  =  7773;  r2  =  0.83).  The  mean  surface 
temperature  for  the  SF  outer  transect  was  1 1.8°C  (<20  m); 
the  mean  temperature  for  the  Segment  2  observations  was 
9.3°C  (<20  m).  In  both  cases,  the  temperatures  were  well 
below  the  recent  and  historical  NDT  calculations.  Although 
strictly  concomitant  nitrate  observations  were  not  obtained 
with  the  SF  fluorescence  and  absorption  data,  it  is  probable 
that  surface  nitrate  was  present  at  detectable  concentrations 
along  both  transects.  Indeed,  the  aforementioned  outlier 


HPLC  and  low-salinity  (S  =  33.28)  Niskin  bottle  sample 
contained  10.6  /iM  nitrate.  Other  macronutrients  (phosphate, 
silicilic  acid)  were  detected  and  generally  abundant  with 
respect  to  nitrate  at  Redfield  et  al  [1963]  proportions  (data 
not  shown).  These  observations  strongly  suggest  that  mac¬ 
ronutrient  limitation  or  stress  is  unlikely  to  explain  the  bio- 
optical  and  fluorescence  trends  in  these  data. 

[32]  The  transition  from  the  pelagic  nanoflagellate  to  the 
neritic  diatom  phytoplankton  communities  in  this  transi¬ 
tional  system  may  instead  be  influenced  by  varying  degrees 
of  iron  limitation  [Ryan  et  al .,  2005;  Bndand  et  al .,  2001; 
Fitzwater  et  al ,  2003;  Johnson  et  al ,  1999;  King  and 
Barbeau ,  2011].  Supporting  evidence  for  this  iron  hypothe¬ 
sis  would  be  a  corresponding  trend  in  surface  dissolved  iron 
concentrations  from  low  values  (<1.0  nM)  at  lower  salinities 
offshore  to  higher  values  (>1.0  nM)  corresponding  to  higher 
salinities  over  the  continental  shelf.  Indeed,  this  pattern  was 
observed  during  the  MUSE  (2000)  Monterey  Bay  field 
experiment  (9-31  August;  Figure  13).  Surface  dissolved  iron 
concentrations  in  excess  of  1  nM  were  observed  at  surface 
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Figure  9.  (a)  Surface  Scalar  PAR  as  measured  by  the  QSR  2200  sensor  on  board  the  R/V Point  Sur  dur¬ 
ing  the  outer  transect,  (b)  The  corresponding  depth  distribution  of  calculated  apparent  variable  fluores¬ 
cence  values  (F*v);  estimated  in  situ  downwelling  planar  PAR  irradiance  values  are  indicated  for  400, 
150,  67,  and  15  fiE  m“2  s_1.  (c)  The  corresponding  surface  (<5  m)  salinity  and  chl-ap  observations. 


salinities  above  33.4,  but  not  at  salinities  below  this  thresh¬ 
old  value  (Figure  12).  Other  iron  data  collected  in  this  region 
confirm  this  trend  [see  Johnson  et  al. ,  2001].  Recall  that  the 
33.4  surface  salinity  isopleth  has  been  previously  used  to 
delineate  subarctic-origin  surface  waters  of  the  California 


Current  from  those  of  the  coastal  transition  zone  [Collins  et  al ., 
2003].  Culture  studies  of  isolated  phytoplankton  strains  from 
the  North  Pacific  and  suggest  iron-limitation  of  the  specific 
growth  rate  for  large  diatoms  (>20  /m\  cell  diameter)  occurs 
at  ambient  dissolved  iron  concentrations  below  1  nM  [ Takeda , 
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Figure  10.  (a)  Apparent  variable  fluorescence  and  (b)  apparent  fluorescence  efficiency  versus  irradiance 
(as  in  Figure  8)  for  a  shelf  transect  (ST)  northwest  of  MB  (12:00-12:53  local  time;  see  Figure  lc).  The 
smoothed  trend  line  (red)  inflection  points  are  indicated  by  the  dashed  lines.  The  blue  line  is  the  trend  from 
the  offshore  transect  (OT;  Figure  8),  for  comparison. 


2011].  In  contrast,  pelagic  phytoplankton  isolates  appear 
better  adapted  to  low-dissolved  iron  conditions  (<1  nM)  con¬ 
ditions  [Takeda,  2011;  Strzepek  and  Harrison ,  2004]. 

4.  Discussion 

[33]  The  transition  from  the  comparatively  low-dissolved 
iron  environment  of  the  subarctic-origin  California  Current  to 
the  more  iron-replctc  upwclling  waters  along  the  California 
coast  is  attendant  to  a  phytoplankton-species  composition 
transition  (pelagic  nanoflagellates  to  diatomaceous  micro¬ 
phytoplankton)  that  is  manifest  as  not  only  a  change  in 
accessory  pigment  composition  (Table  1)  but  also  a  probable 
shift  in  the  architecture  of  the  photosynthetic  apparatus.  It  is 
this  latter  photophysiological  shift  that  may  explain  the 
principal  results  from  sections  3.1  and  3.2:  a  two  to  fourfold 
increase  in  the  apparent  fluorescence  per  unit  absorption 
(AFEp)  or  chi  (AFEC)  and  a  simultaneous  halving  of  the 
apparent  irradiance  saturation  (I'k;  Figure  10b).  For  example, 
these  Tk  (/*E  m~2  s"1)  values  are  directly  related  to  other 


photophysiological  phytoplankton  variables  [Kolber  and 
Falkowski ,  1 993;  Kiefer  and  Reynolds,  1992] 


where  <7PSii  is  the  effective  (functional)  absorption  cross- 
section  of  PSII  (m2  p E_l  for  unit  consistency  but  more 
commonly  expressed  as  A2  quanta'  *)  and  rP  (seconds)  is  the 
minimum  steady  state  transit  time  for  an  electron  to  transition 
from  water  to  the  terminal  electron  acceptor  at  light  satura¬ 
tion.  The  effective  absorption  cross-section  is  the  product  of 
the  optical  absorption  due  to  photosynthetic  pigments  and 
the  probability  that  an  cxciton  will  be  transferred  to  the 
PSII  reaction  center  core  to  perform  photochemical  work 
[Mauzerall  and  Greenbaum ,  1989].  Absorption  cross-sections 
arc  more  generally  a  measure  of  the  effective  target  area  for 
photons  to  strike  the  photosynthetic  unit.  Larger  PSII  “target 
areas”  may  manifest  as  a  higher  apparent  yield  of  fluoresce  per 
unit  chi  since  ~95%  of  fluorescence  per  unit  chi  signal 
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Table  1.  Physical  and  HPLC  Data3 


Mean  Standard  Deviation  Value  Z-Scorc 


Salinity 

33.92 

0.10 

33.28 

-6.40 

Density 

25.96 

0.15 

25.29 

-4.47 

[But  fuco] 

0.010 

0.01 1 

0.074 

5.81 

[Hex  fuco] 

0.034 

0.028 

0.217 

6.48 

[Lutein] 

0.003 

0.002 

0.016 

5.42 

[PPC]/[TCaro] 

0.776 

0.064 

0.590 

-2.90 

[PSC]/[TCaro] 

0.224 

0.064 

0.410 

2.90 

[PPC]/[TPig] 

0.075 

0.024 

0.160 

3.50 

[PSC]/[TPig] 

0.925 

0.024 

0.840 

-3.50 

Microphytoplankton  Fraction 

0.904 

0.091 

0.300 

-6.62 

Nanophytoplankton  Fraction 

0.057 

0.064 

0.490 

6.81 

Picophytoplankton  Fraction 

0.040 

0.031 

0.210 

5.57 

[TChl  b]/[TChl  a] 

0.016 

0.013 

0.094 

6.04 

[Hex  fuco]/[Fuco] 

0.031 

0.142 

1.080 

7.41 

“The  low-salinity  surface 

sample  is 

compared  to  the  total 

sample 

distributions  (n  58)  for  HPLC  variables  where  the  S  <  33.3  sample  meets 
the  p  <  0.01  (|z|  >  2.58)  criterion.  Note:  19'-butaloyloxyfucoxanthin  = 
[But  fuco];  19'-hcxanoyloxyfucoxanthin  =  [Hex  fuco];  PPC  -  photoprotcctivc 
carotenoids;  PSC  =  photosynthetic  carotenoids;  TCaro  *  total  carotenoids; 
TChl  b  =  total  chlorophyll  b  pigments;  TChl  a;  total  chlorophyll  a  pigments; 
Fuco  *  fucoxanthin. 

emanates  from  PS  II  [Johnson  et  al ,  1997;  Kiefer  and 
Reynolds ,  1992],  whereas  ehl  is  distributed  within  PSII  and 
photosystem  1  (PS1). 

[34]  Clearly  a  doubling  of  <tPSh  would  also  likely  result  in 
a  halving  of  1^  via  equation  (5).  Change  in  o>su  by  a  fac¬ 
tor  of  ~1.5  to  3  is  within  the  range  of  natural  variability 
for  crpsn  [Suggett  et  al .,  2004,  2009].  Active  fluorescence 
techniques  (FRRF)  are  routinely  used  to  estimate  rTPSii  via 
examination  of  the  increase  in  fluorescence  (F)  from  mini¬ 
mal  (F0)  to  maximal  (Fm)  under  successive  flashes  of  light. 
For  example,  Moore  et  al.  [2006]  found  a  linear  correspon¬ 
dence  between  the  ratio  of  the  photosynthetie  carotenoid 
1 9T>ut  to  chlorophyll-tf  (HPLC  data)  and  rrpsj,  (FRRF)  where 
surface  crPS1,  values  appeared  to  double  following  a  hori¬ 
zontal  cross-shelf  gradient.  Similarly,  the  outlier  low-salinity 
HPLC  sample  in  these  data  (Table  1)  contained  a  19T)ut: 
ehlorophyll-0  pigment  ratio  an  order  of  magnitude  above 
(0.072  g:g)  the  higher-salinity  19'but:chlorophyll-a  average 
(0.003  g:g).  Thus  the  difference  in  taxonomic  composition 
implied  by  this  particular  HPLC  contrast  (pclagophytcs) 
may  as  well  be  indicative  of  taxon-speeifie  changes  in  <tPSu. 

[35]  Further  extrapolation  from  our  HPLC  data,  however, 
is  problematic  since  the  sample  was  taken  from  the  surface 
(<5  m  depth)  and  contained  comparatively  elevated  photo- 
protective  carotenoids  (Table  1).  The  significant  decline  in 
(negative)  F*v  for  low-salinity  surface  waters  observed 
elsewhere  (Figure  8),  suggests  severe  light  stress  and  so 
photoprotcctivc  carotenoids  in  the  light  harvesting  com¬ 
plexes  for  this  sample  may  have  served  to  enhance  the 
thermal  dissipation  pathway  at  the  expense  of  lowering 
flpsn-  Conversely,  the  generally  lower  rrPsn  for  large  phy¬ 
toplankton  may  serve  as  a  competitive  advantage  in  a  high 
light  and  resource  rich  environment  nearer  the  coast:  smal¬ 
ler  antenna  complexes  may  be  more  rapidly  repaired  and 
reconstituted  from  damage  caused  by  excessive  irradiancc 
[Falkowski  and  Oliver ,  2007].  Indeed,  shelf  observations  of 
F*v  at  800  fjE  m  2  s_1  are  approximately  the  same  as  at 
10  fiE  m~2  s-1  (Figure  10);  this  is  at  least  suggestive  of  a 
superior  acclimativc  capacity  under  the  conditions  of  super¬ 
saturating  irradiance. 


[36]  The  other  term  in  equation  (5),  the  photochemical 
electron  transport  rate  (l/rP),  is  subject  to  similar  variations 
as  a  result  of  light  acclimation/stress  as  well  as  macro-  or 
micronutricnt  deficiencies  [Falkowski  et  al.,  1992].  How¬ 
ever,  recent  work  has  highlighted  some  broadly  applicable 
taxon-specific  trends  in  oPSu  that  may  explain  these  data. 
Suggett  et  al.  [2009]  provide  a  comprehensive  review  of 
laboratory  and  field  FRRF  <rP Su  estimates;  they  find  that  large 
changes  in  aPsll  are  most  significantly  correlated  to  taxon- 
speeific  gradients  in  cell  size  (larger  cells  — ►  smaller  <tPSh). 
This  trend  was  largely  confined  to  eukaryotes,  but  it  is  cer¬ 
tainly  consistent  with  the  nanoflagellatc  to  large  diatom  trend 


Chi  -a  mg  nr3 

Figure  11.  (a)  Historical  data  (as  referenced  in  the  main 
text)  for  phytoplankton  cell  counts.  Epi fluorescence  stained 
cell  volumes  were  converted  to  carbon  biomass  estimates 
using  literature  values,  (b)  These  carbon  estimates  were  con¬ 
verted  here  to  chlorophyll  estimates  using  a  mean  carbon-to- 
chlorophyll-a  ratio  of  47.2.  The  photoautotrophic  taxonomic 
data  were  grouped  into  diatoms,  picophytoplankton,  and 
then  all  others  (by  difference  of  diatoms  and  pieophyto- 
plankton  from  the  total  photoautotrophic  biomass).  The 
resulting  percentages  of  the  total  were  sorted  by  total  chloro¬ 
phyll  concentration  and  smoothed  by  a  moving  10-point 
mean. 
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Figure  12.  (a)  Historical  data  linear  regression  for  nitrate 
versus  temperature  measurements.  Observations  were 
restricted  to  the  upper  200  m  and  latitude  east  of  —123  W. 
(b)  The  same  nitrate  versus  temperature  regression  for  the 
BIODSACE  2008  field  survey. 

identified  in  this  MB  system.  Some  of  this  cell  size  trend  may 
be  due  to  pigment  packaging:  larger  cells  “package”  pig¬ 
ments  more  densely  thereby  reducing  the  effective  target  area 
for  photons  to  strike.  It  is  also  true  that  pigment  packaging 
may  have  contributed  to  the  reduction  in  the  AFEC  values 
observed  nearshore:  the  near-red  chlorophyll-a  absorption 
band  overlaps  with  the  fluorescence  emission  band  and  may 
contribute  to  fluorescence  re-absorption  within  the  cell  [Huot 
and  Babin ,  2010]. 

[37]  Nevertheless,  two  compelling  pieces  of  evidence 
strongly  suggest  that  this  package  effect  may  be  secondary 
to  a  more  broadly  applicable  taxonomie  shift  in  <TpSn.  First, 
different  values  of  <7PSn  for  a  mixed  phytoplankton  popula¬ 
tion  tend  to  mix  conservatively  and  may  be  predieted  using 
simple  mixing  ratios  based  on  the  fluorescence  properties  of 


the  predominant  taxa  [Suggett  et  al ,  2004].  This  property  of 
fluorescence  properties  in  mixed  phytoplankton  populations 
may  explain  the  “quasi-eonservative”  mixing  of  AFEc  evi¬ 
dent  when  plotted  as  a  function  of  surface  salinity  (Figure  6), 
i.e.,  it  is  simply  a  consequence  of  the  mixing  of  pelagic 
nano  flagellates  (high  AFEc,  low  lie,  and  potentially  large 
0>sn)  lafge  coastal  diatoms  (low  AFEc,  high  ITc,  and 
potentially  smaller  <jpsn).  Second  and  as  alluded  to  previ¬ 
ously,  prevailing  depletion  of  essential  micronutrients  may 
place  selective  evolutionary  pressure  toward  larger  <7PSn  in 
pelagic  species.  It  is  far  less  iron-intensive  to  focus  light 
harvesting  capability  upon  expansion  of  antenna  pigment 
complexes  (requiring  nitrogen)  than  it  is  to  replicate  addi¬ 
tional  photosynthetic  reaction  centers  given  their  compara¬ 
tively  expensive  iron  requirements  [Strzepek  and  Harrison , 
2004].  Given  that  the  source-water  region  for  the  California 
Current  is  one  of  the  major  High  Nutrient  Low  Chlorophyll 
(HNLC)  zones  in  the  world  ocean  [Martin  et  al .,  1989], 
it  would  not  be  surprising  that  the  flora  endemic  to  these  sur¬ 
face  waters  should  exhibit  a  comparatively  iron-conserving 
photosynthetie  architecture. 


5.  Conclusions 

[3«]  Consistent  with  the  goal  of  the  BIOSPACE  pro¬ 
gram,  this  analysis  indicates  that  the  fundamental  bio- 
optieal  and  ecological  distinction  between  coastal  MB 
and  California  Current  surface  waters  is  the  transition  from 
mierophytoplankton  eutrophy  to  nanoflagellate  mesotrophy. 
A  coupled  ecosystem-hydrodynamic  model  of  the  region 
would  need  to  resolve  this  transition  in  a  way  that  maintains 
fidelity  to  the  observed  physical  and  bio-optical  property 
features  of  speeifie  water  masses.  This  analysis  also  suggests 
that  there  exists  a  very  clear  photophysiological  distinction 
between  these  respective  phytoplankton  communities  that 
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Figure  13.  Surface  dissolved  iron  data  versus  ship  under¬ 
way  surface  salinity  measurements  (R/V  New  Horizon ) 
obtained  during  the  MUSE  experiment  in  the  Monterey 
Bay  region  ( 1 9-3 1  August  2000). 
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may  be  expressed  as  a  marked  difference  in  irradiancc  satu¬ 
ration  values  and  potentially  the  effective  absorption  cross- 
section  of  PSII.  Circumstantial  evidence  from  literature  and 
historical  data  indicate  that  prevailing  micronutricnt  condi¬ 
tions  may  be  a  contributing  factor  to  this  distinction;  further 
studies  involving  collocated  in  situ  variable  fluorescence 
measurements  (FRRF)  and  dissolved  iron  measurements  are 
required  to  determine  these  relationships  more  conclusively. 


Appendix  A 

Al.  Calculation  of  In-Water  Downwelling  Irradiance 

[39]  The  calculation  of  downwelling  planar  irradiance 
(Ed(A))  at  1  nm  spectral  intervals  was  performed  for  each 
profile  and  each  depth  increment  via  an  iterative  hyperspec- 
tral  attenuation  scheme.  Incident  PAR  was  taken  from  the 
nearest  space/timc  ship  recording  of  surface  PAR  via  the 
R/V Point  Sur  Underway  Data  Acquisition  System  (UDAS); 
specifically,  the  Biospherical  QSR  2200  SPAR  sensor.  This 
is  a  broadband  sensor  and  so  the  spectral  decomposition  over 
the  400-700  nm  range  was  accomplished  via  application  of 
a  spectral  atmospheric  irradiance  transmission  model  [Gregg 
and  Carder ,  1990],  Hyperspectral  percentages  of  the  total 
PAR  output  were  applied  to  the  QSR  2200  data  for  a  full 
spectral  decomposition  estimate  for  the  incident  irradiance 
values.  The  scalar  QSR  2200  irradiance  data  were  converted 
to  downelling  planar  irradiance  by  a  constant  factor  of  0.8 
[Beardsley  and  Zaneveld,  1969;  Jerlov ,  1976].  This  was  an 
approximation  and  the  ratio  Ed/E0  varies  spectrally  and  as  a 
function  of  the  angular  distribution  of  the  light  field.  Broad¬ 
band  air/sca  boundary  reflectance  was  held  constant  at  4%. 

[40]  Hyperspectral  downwelling  planar  irradiance  (Ed(A,  z)) 
was  then  attenuated  iteratively  through  the  water  column 
at  (z)  depths  following  the  single  scattering  approximation 
[Sathyendranath  and  Platt ,  1 989]  for  the  diffuse  attenuation 
at  each  depth  and  wave  band 


Kd(A,z)  = 


a{\z)  +  bb(\,z) 

PdW 


(Al) 


The  spectral  range  for  the  mean  cosine  of  downwelling 
irradiance  (//d(A))  was  estimated  from  published  tabular 
values  [Morel  and  Maritorena ,  2001],  and  confined  to  a 
range  of  0.767-0.901.  Total  absorption  coefficients  from  the 
AC-9  instrument  at  8  channels  [412,  440,  488,  510,  532, 
555,  650, 676  nm]  and  spectral  backscattcring  coefficients  at 
3  channels  [470,  532,  660  nm]  from  the  ECO-BB3  sensor 
were  used  to  interpolate  hyperspectral  diffuse  attenuation 
values  at  1  nm  increments.  Attenuation  through  the  water 
column  of  a  spectrally  decomposed  diffuse  light  field  may 
be  approximated  by  the  Beer-Lambert  Law 


Ed(A,z)  =  Ed(A,z-l)e(-MA-')A').  (A2) 


[41]  Iterative  application  of  A2  was  used  to  attenuate  each 
wave  band  ( 1  nm)  through  the  water  column.  Irradiance  at 
each  depth  was  then  calculated  as  the  sum  of  the  spectrally 
decomposed  irradiance  over  the  PAR  spectral  range  (400- 
700  nm). 


A2.  MUSE  Dissolved  Iron  Data  and  MBARI 
Mooring  Data 

[42]  Data  on  surface  dissolvable  iron  were  obtained  from  the 
MOOS  Upper-Water  column  Science  Experiment  (MUSE) 
(19-31  August  2000)  web  site  (http^/www.mbari.org/MUSE/). 
Surface  dissolvable  iron  determinations  were  converted  to 
dissolved  iron  following  methods  described  in  Johnson  et  al 
[2001].  The  processed  dissolved  iron  data  were  matched  by 
time  and  location  to  the  underway  CTD  data  (UCTD  data) 
obtained  from  a  Seabird  SBE-21  (SN  2960)  Conductivity 
Temperature  sensor  mounted  on  the  R/V  New  Horton .  The 
resulting  matchup  between  dissolved  iron  determinations  and 
underway  salinities  are  plotted  in  Figure  13.  Spatial  maps  of 
the  dissolved  iron  distribution  arc  available  at  the  MBARI 
MUSE  web  site.  Historical  and  recent  temperature  and  salinity 
data  from  MBARI  permanent  moorings  are  also  available 
from  the  MBARI  web  site  (www.mbari.org/oasis/). 
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